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TRANSONIC SINGLE-MODE FLUTTER AND BUFFET OF A LOW 


ASPECT RATIO WING HAVING A SUBSONIC AIRFOIL SHAPE 

Larry L. Erickson 
Ames Research Center 


SUMMARY 


Transonic flutter and buffet results obtained from wind-tunnel tests of a low aspect ratio 
semispace wing model are presented. The tests were conducted to investigate potential transonic 
aeroelastic problems of vehicles having subsonic airfoil sections (e.g., the space shuttle). The model 
employed NACA 00XX-64 airfoil sections in the streamwise direction and had a 14° leading edge 
sweep angle. Aspect ratio, taper ratio, and average thickness were 4.0, 0.35, and 8 percent, respec- 
tively. The model was tested at Mach numbers from 0.6 to 0.95 at angles of attack from 0° to 15°. 

Two zero lift flutter conditions were found that involved essentially single normal mode 
vibrations. With boundary layer trips on the model, flutter occurred in a narrow Mach number range 
centered at about Mach 0.90. The frequency and motion of this flutter were like that of the first 
normal mode vibration. With the trips removed flutter occurred at a slightly higher Mach number 
but in a mode strongly resembling that of the second normal mode. 

Several instances of low torsional damping were observed at high angles of attack, but stall 
flutter was not encountered, except possibly at M = 0.90, a » 13°. Maximum buffet loads occurred 
at angles of attack between 10° and 14°. Quantitative results for the maximum wing root bending 
moments due to buffet are presented in a form suitable for scaling. 


INTRODUCTION 


Wing and tail surfaces considered for the space shuttle orbiter have relatively thick airfoil 
sections that are generally associated with subsonic flight (ref. 1 ). Wings having such airfoil sections 
are not normally flown above M cr \ as a result little information is available concerning the transonic 
flutter behavior of thick wings. What information there is has been obtained from wind-tunnel tests 
of aeroelastic models since the unsteady viscous, mixed-flow problem associated with transonic 
flutter is not presently solvable by analytic means. 

Flutter and buffet results from an aeroelastic model of one of the early Manned Spacecraft 
Center straight wing concepts are reported in reference 2. This wing had NACA 00XX-64 airfoil 
sections, an average spanwise thickness ratio of 1 2 percent, and an aspect ratio of 7. Wind-tunnel 
tests of the model revealed an unusual zero lift flutter behavior that occurred over a narrow 
transonic Mach number range, at quite low dynamic pressures, and at a frequency nearly identical 
to the first mode still-air vibration frequency. The present investigation was undertaken to deter- 
mine if a thinner, lower aspect ratio wing model having the same airfoil shape would experience the 


same type of dynamic instability. In addition, tests were conducted at sufficiently large angles of 
attack to determine the stall flutter and buffet characteristics of the model. 


WIND TUNNEL 


The tests were conducted in the Ames 11- by 11 -Foot Transonic Wind Tunnel. This facil- 
ity is a closed return slotted-throat, variable-pressure tunnel with an operating range of 
51 220 kN/m 2 (1 5 <p t < 65 in.-Hg), where p t is the total (stagnation) pressure. The 

variations in free-stream dynamic pressure and Reynolds number with tunnel total pressure are 
shown in figures 1 and 2, respectively. Pt> in _ Hg 
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Figure 1.- Variation in free-stream dynamic pressure with total pressure. 
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Figure 2.— Variation in free-stream Reynolds number with total pressure. 
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GEOMETRY AND CONSTRUCTION 


Figure 3 shows the wing planform geometry. Wing construction details are illustrated in fig- 
ure 4. The airfoil sections were NACA 00XX-64 profiles parallel to the flow with maximum thick- 
ness ratios varying from 9-1/3 percent at the root to 6-2/3 percent at the tip, giving an average 
thickness of 8 percent. The load-carrying structure consisted of a fiberglass-epoxy skin laid over a 
lightweight foam interior, and a birch spar at the 20 percent chord position. The spar and a skin 
thickness change at 30 percent chord were incorporated in the wing to adjust the chordwise posi- 
tions of the section elastic axis and center of gravity. Spar dimensions and skin thicknesses varied 
linearly along the semispan between the values given at y/L = 0 and 1 . 



Figure 3 - Wing planform geometry. 


Figure 4.- Airfoil section construction. 


The model was mounted on a tapered steel tongue that matched a tapered cavity in the base of 
the model. A schematic of the wing support system is shown in figure 5. When the tongue was fully 

inserted into the cavity it fit snugly against the mating r \ -» 

surfaces except at six slots, 2.54 cm (1 in.) wide. These slots / \ 

provided clearance for wire-resistance strain-gage bridges that / \ ... Lending edge 

were bonded to the tongue at y= y/L = 0 .\ 8 . Separate / [ To ~~ o~~o~ T tT \ 

bridges were aligned for both torsion and bending measure- / ! | „ ri n V'A 

ments. Figure 6 shows the model mounted in the wind / j | | j | j ] j \ \ \ 
tunnel. The large fiberglass fairings at the wing root covered I l' ° i !° I !°! I ° '\ \ \ 

clamping bolts used for tightly securing the wing to the 

tongue. ' Flats for Clomping bolts 


Metallic potting maten 


Except where noted, all testing was performed with 
boundary-layer trips located 5.7 cm (2-1/4 in.) from the 

J J , / / 'Cavity for tongue support 

model leading edge. The trips were 3.2 mm (1/8 in.) wide -'slots for strain gQg es 

strips of 0.25 mm (0.010 in.) diameter glass spheres and there 190965 b0 " aed ,0 ,0 " 9 “ el 

were approximately 40 spheres per cm (100/in.) of trip 

length. The model was covered with a white mixture of Figure 5.- Schematic of wing support 
highly viscous oil and titanium dioxide to indicate flow pat- system, 

terns. Flow patterns observed at a = 0°, M = 0.6, showed that the flow was nearly two-dimensional 
over the wing surface and was apparently little affected by the bolt cover fairings. 


Figure 5.- Schematic of wing support 
system. 
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(b) 3/4 rear view. 

Figure 6.— Model mounted in Ames 11- by 11 -Foot Transonic Wind Tunnel. 



Elastic and Mass Properties 


The model was designed to have the same scaled strength as the AR = 7 model described in 
reference 2. The spanwise variation of the section bending stiffness (El) and the bending-to-torsion 
stiffness ratio (EI/GK) for the model is shown in figure 7 for sections parallel to the flow. These 
values for El and GK were computed from the 
wing cross section geometry and using 
E= 14 GN/m 2 (2.0 X 10 6 lbf/in. 2 ) for the fiber- 
glass skin and birch spar and G = 4.1 GN/m 2 
(0.60 X 10 6 lbf/in. 2 ) for the fiberglass. The rigid- 
ity of the foam and the torsional rigidity of the 
spar were neglected. Using these stiffnesses, span- 
wise deflections due to a tip load, and spanwise 
angles of twist due to a tip-torque were computed 
from engineering beam theory. 1 These bending 
and twisting displacements are shown in fig- 
ure 8(a) and (b), respectively. Measured deflec- 
tions obtained from static bending and torque 
tests are also shown. Angles of twist were mea- 
sured in planes perpendicular to the 40 percent 
chord line. The bending deflections were mea- 
sured at 40 percent chord positions with the tip load at 40 percent chord. The fairly good agree- 
ment between the measured and computed deflections indicate that the model stiffness properties 
given in figure 8 are reasonably accurate. 



Figure 7.— Spanwise variation of elastic stiffnesses. 


Additional model section properties, computed from the chordwise distribution of mass and 
bending stiffness, are given in table 1 ( c is the local chord length). The mass ratio p of the portion of 
the model outboard of the tongue is 


16 /310 X 10~ 4 
Poo V Poo 


where p ^ is the free-stream air density in 
kg/m 3 (slugs/ft 3 ). The mass ratio is plotted 
versus Mach number in figure 9 for 
= 51 kN/m 2 (15 in.-Hg). Also shown is the 
variation in first mode reduced frequency, 
k x = bu> 1 /V 00 . (Model values of pand k are 
required to predict full-scale aerodynamic 
damping as explained in the section on 
buffeting.) 


TABLE 1.- MODEL SECTION PROPERTIES 


I 

y = 0.255 

y = 1.0 

Distance from leading edge 

mm 

Ml 

to elastic axis 



to center of gravity 

Ha 

Ha 

Radius of gyration 



about elastic axis 

0.290c 

0.292c 


1 Due to the low aspect ratio and plate like geometry of the model, the bending deflections were computed by 

replacing El with £7/(1 — n 2 ), v = 0.3. Also, equation (g) of reference 3 (p. 65) was used to correct approximately 

the twisting deflections for wing root warping restraint. (This correction gave a 15 percent decrease in angle of twist 

at the wing tip.) 
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Still-Air Vibration Characteristics 


Values of still-air frequencies and damping ratios are given in table 2 for the first three normal 
vibration modes of the model. Mode 1 , a beam-like fundamental bending mode, is shown in 
figure 10(a) being driven at reasonance by an air-jet exciter. Modes 2 and 3 are plate-like modes as 
shown by the node line patterns in figures 10(b) and 10(c), respectively. 


TABLE 2.- FREQUENCY AND DAMPING RATIO 



Strongback installation 

Wind tunnel 
installation 

Mode number 

Frequency, 

Hz 

Damping 

ratio, 

Frequency, 

Hz 

Damping 

ratio, 



percent 


percent 

1 

72.5 

0.5 

70.2 

0.75 

2 

160 

0.95 

— 

— 

3 

205 

0.88 


— 




(b) Second mode;/= 160 Hz. (c) Third mode;/= 205 Hz. 


Figure 10.— Concluded. 

TEST PROCEDURE AND DATA REDUCTION 

Previous tests (ref. 2) of the AR = 1, t/c = 12 percent flutter model having the same airfoil 
profile as the present one revealed a very unusual type of transonic flutter, which occurred only 
over a narrow Mach number range. Because of this flutter sensitivity to Mach number, the present 
test was conducted in the following manner: 

1 . Zero-lift flutter search.— With a = 0°, and the tunnel at the lowest operating pressure, Mach 
number was varied between 0.6 and 0.95. (Limit amplitude flutter was encountered in the vicinity 
of M = 0.90, but the model had sufficient strength to allow the flutter region to be traversed 
without damage.) Total pressure was then increased and the Mach number survey repeated. 

2. Stall flutter and buffet test (a> 0°).— For a> 0°, all tests were conducted at the lowest 
tunnel operating pressure so as to keep static airloads to a minimum. At fixed Mach numbers, the 
angle of attack was varied from 0° to 15°. 

The above procedure was followed with boundary-layer trips in place. After completion of the 
stall flutter and buffet tests these trips were removed and the zero lift flutter search was repeated at 
a low tunnel pressure. With the trips removed, no evidence of flutter was found near M = 0.90 so 
pressure was then increased while Mach number was held between 0.90 and 0.93. Eventually, a 
pressure was reached at which flutter was again encountered. As discussed later, this flutter was of a 
different nature than that found with the trips in place and the model was destroyed. 

The wing motion was continuously monitored by oscilloscope displays of the strain-gage 
signals and by a television display from a camera mounted above the wing tip. High-speed motion 
pictures of the wing motion were taken at several test conditions. 

Throughout the test, the strain-gage signals were continuously recorded on magnetic tape. 
During fixed test conditions, the time histories of the model response were recorded over a 135-sec 
period. The dc components from the torsion and forward bending gage were filtered electronically 
before being recorded; this made the entire range of the recording system usable for the model 
dynamic response. The complete signal (ac + dc) from the rear bending gage was recorded so that 
comparisons could be made of the dynamic and static bending loads on the model. 
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The recorded data were analyzed from measurements of the overall rms, frequency spectra, 
autocorrelations, and oscillograph records of the time histories. 

The still-air frequency and damping ratio of the first bending mode were measured several 
times during the test and were unchanged. Wind-on damping ratios were obtained from envelopes of 
the autocorrelation functions (ref. 4). 


RESULTS AND DISCUSSION 


Zero-Lift Flutter 

At a = 0° the model was tested at total pressures ranging from 51 kN/m 2 to 186 kN/m 2 
(15 in.-Hg to 55 in.-Hg) at Mach numbers from 0.6 to 0.95. Over this pressure range, the model was 
flutter free except in a narrow Mach number range in the vicinity of M = 0.90. The characteristics of 
this flutter are described below. Except where noted, all results are for the model with boundary- 
layer trips attached. 

Response levels — Figure 1 1 shows the effect of Mach number on the rms signal level og of the 
bending gage at the lowest tunnel operating pressure. Evidence of flutter at a = 0° is indicated by 
the large rise in response near M = 0.90. 

Note that this response is completely 
eliminated with the model at a = 2° ? At 
higher pressures, the general character of 
the response was the same as shown in 
figure 12, but the magnitude of og at 
M = 0.90 decreased as p t increased. The 
same behavior was observed on the thicker, 
higher aspect ratio model of reference 2. 


Frequency spectrums.— Frequency 
spectrums of the bending and torsion gage 
signals are shown in figure 1 2 for a = 0° 
and 2° at M - 0.90. At a = 0°, both the 
bending and torsion gages show a pre- 
dominant frequency component at about 
72 Hz, a value nearly identical to the still-air first mode vibration frequency of the model. At 
a. = 2°, the 72 Hz component still predominates in the bending-gage spectrum (but at a much lower 
level), and it is barely noticeable in the torsion-gage spectrum. 

Damping.— Net damping levels (aerodynamic plus structural) associated with the bending 
motion of the first mode frequency are shown in figure 13 where $\g, the ratio of damping to 
critical damping, is plotted against Mach number (p ^ - 51 kN/m 2 (15 in.-Hg)). These damping ratios 
were obtained by passing the bending gage signals through a 65- to 85-Hz bandpass filter and 
computing the autocorrelation of the resulting signal. For a = 0°, the net damping is seen to be less 


2 A careful Mach number survey at a = 2° (0.85 <, M < 0.95) showed no indication of flutter. 
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Frequency spectra of strain gage signal, M d / Hz 


Bending gage 


Torsion gage 




10° r 


a = 2° 


, 72 Hz 


6 | J 

0 250 500 

Frequency, Hz 


Figure 12.— Frequency content of strain gage signals; 
p t = 51 kN/m 2 (15 in.-Hg); = 0.90. 



Figure 13.— Effect of Mach number and angle of attack on 
damping; p f = 51 kN/m 2 (15 in.-Hg). 


Note' amplitude is not drawn to scale. 


= .90 



t = 0 t = P/4 t = P/2 t = 3 P/4 t = P 


Figure 14.— Flutter motion with boundary -layer 
trips-attached. 


Elastic axis 



Center of gravity 


7777777 


Sketch (a) 


than the still-air value in the narrow Mach 
number range centered at about M = 0.90. 
Thus, the large response level of figure 1 1 is 
due to negative aerodynamic damping 
(i.e., flutter) rather than to a severe buffet 
type input. Ranges in damping for 
2 0 <a<6° are also shown in figure 13. 
For a = 2° the aerodynamic damping 
decreases in the vicinity of M = 0.90 but 
flutter does not occur. 

Wing motion.— The limit amplitude 
flutter motion was established with the aid 
of high-speed motion pictures. One cycle of 
this motion at quarter-cycle intervals is 
illustrated in figure 14. The wing bends and 
twists in phase, the nose-up attitude occur- 
ring when the wing is bent down. (The 
same flutter motion was observed on the 
thicker, higher aspect ratio wing of ref. 2.) 
Measurements from high-speed motion 
pictures of the motion indicate the maxi- 
mum double amplitude wing tip displace- 
ment (measured at 0.4 c ( ) and rotation 
were approximately 2.54 cm (1 in.) and 2°, 
respectively. 

It is of interest to compare the flutter 
motion with the free-vibration motion of 
the two-degree-of-freedom system sketched 
below. If the center of gravity and the elas- 
tic axis coincide, one of the two normal 
modes of this system vibrates in pure trans- 
lation at a frequency oo^ while the other 
vibrates in pure rotation at a frequency cog. 
However, if the center of gravity is aft of 
the elastic axis (as in the case of the wind 
tunnel model, and as shown in sketch (a)) 
the normal modes contain contributions 
from both degrees of freedom. In this case, 
and for co^ < c oq, the fundamental (lowest 
frequency) normal mode has the same phas- 
ing (between the translational and rota- 
tional degrees of freedom) as the flutter 
motion shown in figure 14 (ref. 5). This 
similarity in phase, together with the nearly 
identical values of the flutter and funda- 
mental normal mode frequencies of the 
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model, strongly suggests that the motion-induced aerodynamic forces are able to feed energy into a 
single normal mode. (This is in contrast to “classical” flutter, which usually involves the participa- 
tion of two or more normal modes.) 

Speculation concerning the flutter mechanism - Several types of single degree of freedom 
flutter have been predicted on the basis of unsteady potential flow theory (ref. 6). Other types, 
however (e.g., control surface buzz (ref. 7), stall flutter (ref. 8), and hammerhead configuration 
instabilities (ref. 9)), are generally attributed to nonpotential flow characteristics such as oscillating 
shock-wave boundary layer interaction and unsteady separated flow. 

The single mode-like flutter behavior in the present investigation is completely analogous to 
that of the 12 percent thick, aspect ratio = 7 model described in reference 2. The instability of the 
12 percent thick model, which occurred at 0.84 <M < 0.86, appeared to be associated with the 
onset of shock stall. (Oil flow photos indicated a large region of shock-induced separated flow at 
M 0.85, the extent of which was extremely sensitive to small (0.01) changes in Mach number.) 
This possibility is consistent with the present results since the 8 percent thick model would be 
expected to encounter shock stall at a higher Mach number. 

The restriction of the instabilities to narrow Mach number ranges is also an indication of shock 
stall onset. At Mach numbers greater than critical, static lift and pitching moment coefficients of 
airfoils similar to those used on the present model (and on the 12 percent thick model of ref. 2) can 
display large variations in magnitude over narrow Mach number ranges (refs. 10 and 11). The 
chordwise static pressure data of reference 1 1 indicate these variations are most extreme at the 
Mach number corresponding to the onset of shock -induced trailing edge separation (M se p). At this 
condition (i.e., where the flow is transitional between attached and separated) there is likely to be a 
strong interdependence between wing motion, shock position and strength, and boundary layer 
separation. For example, the upper and lower surface shocks, in conjunction with the rotational 
motion 3 of the wing (due to a normal mode vibration), might cause the boundary layers on 
opposite sides of the wing to alternate between attached and separated flow. Thus, self-induced 
oscillating forces would be produced in synchronism with the wing motion. With the proper magni- 
tude and phase, these oscillating forces could do enough positive work per cycle to balance the 
energy loss due to structural damping. An instability produced by such a mechanism would be 
expected to occur only in a narrow Mach number range centered about M sep since at slightly lower 
Mach numbers the flow would be predominately attached, whereas at slightly higher Mach numbers 
the flow would be fully separated between the trailing edge and the upstream shock-induced 
separation point. 

In the flutter mechanism postulated above, the fluctuating forces are considered to act in 
synchronism with the motion of a normal mode vibration. This would explain the flutter frequency 
being nearly identical to a normal mode frequency. (Some difference in frequency would be 
expected due to aerodynamic mass and stiffness effects.) However, the ability of a particular mode 
to become unstable probably requires a reasonable fit between the spatial pattern of the normal 
mode motion and the spatial pattern of the fluctuating surface pressures. This implies that normal 

3 Unpublished test results for an aluminum skin, honeycomb core version of the AR = 7 fiberglass model of 
ref. 2 indicate the rotational degree of freedom is very important. The aluminum model was approximately three 
times stiffer in torsion than the fiberglass model but had approximately the same bending stiffness and first mode fre- 
quency. This torsionally stiffer model was flutter free at all test conditions where the fiberglass model experienced 
flutter. 
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modes other than the fundamental could also become unstable by altering the spatial distribution of 
the separated flow. Evidence that such behavior is possible is given in the following section. 


Natural transition.— At a given tunnel condition (fixed Reynolds number per unit length), 
removal of the boundary-layer trip should change the location of boundary layer transition and the 
boundary layer thickness from that of the tripped model. Changes in tunnel Reynolds number 
would produce additional changes in these properties. The attendant effect on the occurrence of, 
and the spatial distribution of, shock-induced flow separation (refs. 12, 13, 14) and the resulting 
forces, might then alter the wing flutter characteristics. Reynolds number variations, for example, 
are known to affect static pitching moments at supercritical Mach numbers (see ref. 15, figs. 4 
and 6). 


o,D, With boundary layer trip 


With the boundary-layer trips removed, a Mach number sweep at tx = 0° from M = 0.6 to 0.95 
at p t = 84 kN/m 2 (25 in.-Hg) showed no evidence of flutter. The change in wing damping due to 
removal of the boundary -layer trips is shown in figure 15. Damping ratios associated with the 

fundamental mode frequency are given at 
two total pressures for the wing with trips 
(open symbols), and at the higher of the 
two pressures for the wing without trips 
(dark symbols). With the trips in place the 
aerodynamic damping becomes negative at 
both pressures, although confined to a 
narrower Mach number range at the higher 
pressure. At this higher pressure, but with 
the trips removed, there is a small decrease 
in damping near M = 0.9, but flutter does 
not occur. Thus, Reynolds number effects 
apparently play a crucial role in the flutter 
mechanism. 



p. = 84 k N/m 2 /''' 

(25 in- Hq) \ 

V 

o^/o A \ 

f 

w \ 1 


p.= 5l kN/m 2 \1 

(15 in.-H g ) M 

\f 

” Still -air damping at atmospheric pressure \| 



if , 


Figure 15.- 


Effect of total pressure and boundary-layer trip on 
damping; a = 0°. 


Since the model without trips was flutter free at p t - 84 kN/m 2 (25 in.-Hg), Reynolds number 
(as well as dynamic pressure) was increased by increasing pp the Mach number being held between 
0.90 and 0.93. Flutter was encountered at p t = 154 kN/m 2 (45.7 in.-Hg), M = 0.92. No flutter 
occurred at this test condition when the boundary -layer trips were in place; thus, the flutter 
encountered without trips is evidently due to a Reynolds number effect rather than to the increas- 
ing dynamic pressure. The model was destroyed at this condition just after several seconds of 
high-speed motion pictures were taken of the flutter motion. 

Two frames from the motion picture film are reproduced in figure 16(a). These frames show 
the model near its extreme flutter amplitudes. At the wing tip, the amplitude of the leading edge 
motion is greater than that of the trailing edge. Also, the trailing edge has a nodal point a short 
distance inboard of the tip section. Figure 16(b) is a sketch of this flutter motion; note that it is 
very similar to the second still-air normal vibration mode of the model shown in figure 10(b). 

The flutter frequency in this case was 145 Hz, a value about 10 percent less than the second 
mode still-air frequency of 160 Hz. The 160 Hz still-air value is evidently decreased by aerodynamic 
mass and stiffness effects, as evidenced by the data in table 3, which gives frequencies corresponding 
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to the predominant peaks in the strain-gage signal spectrums 
for various flow conditions. For M <0.91, only the first 
mode frequency (72 Hz) is predominant. For M > 0.92, the 
second mode frequency is also apparent; at the low total 
pressure it has a value of 160 Hz (same as the still-air value) 
for the model with trips, and a value of 155 IIz for the 
model without trips. At the higher total pressure these 
values have been reduced by 10 Hz. Thus, the 145 Hz flut- 
ter frequency for the model without trips can be traced to 
the model’s second still-air normal mode frequency. This 
finding, together with the similarity between the second 
normal mode motion and the flutter motion, is again indica- 
tive of a single mode type of instability. 



The above experimental results indicate that different 
normal modes can become unstable due to the presence or 
absence of a boundary-layer trip. This supports the conten- 
tion that the ability of a particular mode to become 
unstable is influenced by the spatial distribution of the 
boundary layer characteristics and the attendant effects on 
the onset of shock stall. 

Implication of results.— The results of the present 
investigation and of reference 2 demonstrate that either 
high or low aspect ratio straight wings with subsonic airfoils 
are susceptible to flutter over narrow transonic Mach 
number ranges. The results obtained apply only to the con- 
figurations tested, but they raise concern for other wing and 
tail geometries since the flutter is associated with single 
mode type instabilities rather than with coupled mode inter- 
actions. For example, the node line pattern of the unstable 
second mode for the AR = 4 model (fig. 10(b)) resembles 
that expected for delta wing configurations. 

Several aspects of these instabilities are of secondary 
importance to more usual flutter behavior. For one thing, 
the instabilities can apparently be prevented by flying at 
moderately small angles of attack. Thus, the problem could 
be avoided for a wing (but not for nonlifting tail surfaces). 
Second, the occurrence of flutter in very narrow Mach 
number ranges and the limit amplitude behavior may make 
it possible (but not necessarily acceptable) to fly through 
such flutter regions with no structural damage. Accum- 
mulated fatigue damage could be a major consideration, 
however. 

These Mach number sensitive results indicate that very 
close attention should be paid to small Mach number 



(a) Frames from highspeed motion 
pictures; flow is from top to 
bottom of picture. 



(b) Sketch of flutter mode shape. 

Figure 16.— Flutter motion with 
natural transition; M ^ = 0.92, 
p t = 154 kN/m 2 (45.7 in.-Hg). 
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TABLE 3 - PREDOMINANT FREQUENCIES IN STRAIN-GAGE SIGNAL SPECTRUMS, Hz; a = O' 



With boundary-layer trips 

Without boundary-layer trips 


p t = 84 kN/m 2 

p t = 152 kN/m 2 

p t ~ 84 kN/m 2 

p t = 154 kN/m 2 


(25 in.-Hg) 

(45 in.-Hg) 

(25 in.-Hg) 

(45.7 

in.-Hg) 


Bending 

Torsion 

Bending 

Torsion 

Bending 

Torsion 

Bending 

Torsion 


gage 

gage 

gage 

gage 

gage 

gage 

gage 

gage 

M = 0.85 

72 

None 

72 

None 

— 


_ _ _ 

_ 

0.89 0.91 

72 

72 

72 

72 

72 

None 


— 

M = 0.92 

72;160 

160 

72 a 

150 a 

72 ;1 55 

155 

145 

145 

M = 0.95 

72;1 60 

160 

72;150 

150 

— 


— 

— 


a Value interpolated from frequencies at lower and higher total pressures. 


variations in transonic flutter tests. Because of the demonstrated Reynolds number effects, con- 
sideration also should be given to simulating full scale boundary layer characteristics as closely as 
possible. 


Buffet and Stall Flutter 


The following results describe the dynamic response of the wing for angles of attack from 0 
to 1 5° at Mach numbers from 0.6 to 0.95. All results were obtained at a total pressure of 51 kN/m 2 


(15 in.-Hg) and with boundary-layer trips attached 



Figure 17.— Relative variation in rms bending and torsion 
gage signals with angle of attack; M x - 0.90, 
p ( = 51 kN/m 2 (15 in.-Hg). 


to the model. 

A measure of the dynamic response of the 
wing is given by the standard deviations of the 
bending and torsion gage signals <jg and oj, 
respectively. An example of the relative varia- 
tion in og and aj with angle of attack is shown 
in figure 17, which is for M = 0.90. The large 
response at a = 0° is due to the first mode flut- 
ter condition discussed earlier and does not 
occur at Mach numbers slightly less than or 
greater than 0.90. Except for this M 0.90 case, 
ag and a j remain relatively constant until some 
angle is reached at which a fairly rapid increase 
in response occurs as a increases. As shown in 
figure 17, ag and Oj peak at relatively large 
angles of attack and then diminish as a increases 
further. 


For each test Mach number, examination of the strain-gage frequency spectra at the angles of 
maximum response (a > 0°) shows that ag (bending response) is almost entirely due to vibration at 
the first normal mode frequency, whereas oj- (torsion response) is primarily due to the second 
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I 

I 

mode vibration. (The spectrum peaks at these vibration frequencies were one to two orders of 
magnitude larger than those at any other frequency.) 

Net damping ratios for the first (predominantly bending) and second (predominantly torsion) 
modes are denoted as ^ and respectively, 4 and are shown in figure 18 for M = 0.90. For 
10° 15° has a constant value o5 

several times larger than the first mode still- 
air damping value. This implies that the 
a « 13° peak in (fig. 17) is due to an 
increase in dynamic load input and can be 
categorized as wing buffeting. The peak in s 03 
oj, however, occurs where aerodynamic 
damping $ 2 T near zero - Hence, the peak ~ .02 
in aj is more closely associated with the 
wing being at, or near, a condition of stall 0 
flutter (ref. 1 6). 

; The total wing response at large angles 0 

I of attack can thus be regarded as a combin- 

ation of buffeting and near stall flutter. Figure 18.- Variation in damping ratio with angle of attack; 

However, in this case the buffet (relatively M^-0.90, p ( - 51 kN/m i 2 (15 in.-Hg). 

| large damping) and stall flutter (low damping) phenomena can be fairly well distinguished since 

j they affect different modes, which are well separated in frequency (i.e., response due to modal 

frequency overlap (ref. 1 7) is probably not significant). For convenience, then, the peaks in bending 
1 response (<7g) will be referred to simply as buffeting. 

i A buffet region can be estimated from figure 1 7 and similar figures for other Mach numbers. 

For a given Mach number, the angle of buffet onset was determined as shown in figure 19. 
Figure 20 shows the approximate Mach number and angle of attack combinations at which the 
buffet boundary and maximum buffet intensity occur. At M = 0.90 the buffet onset angle is 

I 

/Maximum buffet intensity 




Figure 20.— Buffet onset boundary and 
angles of maximum buffet intensity; 
P t — 51 kN/m 2 (15 in.-Hg). 


4 The damping ratios, f for modes 1 and 2 were obtained by bandpassing the strain-gage signals about 

the First and second mode frequencies, respectively, and computing the autocorrelations of the resulting signal. 
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approximately 6°, a value near that at which (?i J g) max occurs, (see fig. 18). 5 This behavior has 
been previously noted in reference 18 where f was measured from flight data. Quantitative 
information concerning the points of maximum buffet intensity is given below in a form useful for 
predicting full-scale buffet loads. 

The total wing root bending stress (static plus dynamic) is equal to the static stress times the 
factor [ 1 + (Mdy n /M stat j c )] where M stat { c and M ( j yn are the static and dynamic wing root bending 
moments, respectively. As mentioned earlier, the dynamic bending gage signal o g ^W^ yn is due 
almost entirely to the first mode vibration. Consequently, full scale values of M^y n /M sta ti c can be 
predicted from model values by the following buffet scaling relation (ref. 19): 


^dyn^ static? full scale _ /(? a + )] full scale 

Wdyd M static? model [k t / (£ a + )] model 


where k x = &co 1 /F 00 (first mode reduced frequency) and and are the aerodynamic and struc- 
tural damping ratios (fraction of critical) associated with the first mode vibration frequency. 6 Use 
of the scaling relation for aerodynamic damping due to bending (ref. 20) 


)f u u sca j e (? a M£i) m odel 


( 2 ) 


and the fact that co 2 °c( 1 /£ 4 ) ( EI/m w ) 0 , allows equation (1 ) to be written as 


R = 


M dvn ^ 

Jinyyi 

_ F^dyn 

/ (1 + 

^ static 

| (EIJZ*q) 

model M static \ 

( (EIJZ 4 q) 


full scale 


(3) 


This form may at times be more useful than equation (1) since is often small compared to 
unity. The moment ratio then scales (approximately) by the ratio of a term that is more 

readily estimated than Aq and for the full-scale wing since it depends on static rather than 

dynamic structural properties. 


Because buffeting is a random type process, can fluctuate considerably from its root- 

mean-square value. In the present investigation, the peak values of the dynamic bending moment 
( Mp ea k-dyn ’ big- 21) experienced by the model were determined from bending gage oscillograph 
records representing 120 sec of the time history (« 8500 cycles 7 at the first mode frequency of the 
wing). Model values of R (eq. (3)) corresponding to ^lp ea k-dyn are §i ven i n figure 21 for the Mach 

5 The peak in at a = 6° for M- 0.90 is considered to be coincidental. Comparison of buffet onset angles 
with angles of (? 27 ^max at other Mach number shows no apparent correlation except that (forVnax a bvays 
occurred at angles less than or equal to the buffet onset angles (e.g., 50 percent less at M = 0.60). 

6 The scale relation given by equation (1) requires the model and full scale wings to have similar geometry, and 
similar generalized modal characteristics and reduced frequencies pertinent to the first vibration mode. 

7 For a stationary narrow-band Gaussian process with zero mean, the expectancy of exceeding a j 3cr 1 value 
would be once every 45 cycles (ref. 21). 
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Figure 21.- Peak wing-root bending moment and first 
mode damping at angles of maximum buffet intensity. 



Figure 22.- Variation in second mode torsional damping 
ratio with angle of attack; p f = 51 kN/m 2 (15 in.-Hg). 


number angle of attack combinations of max- 
imum buffet intensity. The first mode damp- 
ing ratios used to compute R are also given in 
this figure (f s was taken as the still-air damp- 
ing value). The model values for $ a can be 
scaled to full scale values by use of equa- 
tion (2) and the value for p and k x given in 
figure 9. 

As mentioned earlier, tendencies toward 
stall flutter instabilities can be distinguished 
from the forced buffet response by examining 
the variation in j with angle of attack. This 
variation was shown in figure 18 for M = 0.9 
and is shown in figure 22 for three additional 
Mach numbers. In each case values of j less 
than the still-air damping value would cor- 
respond to negative aerodynamic damping 
and a condition of stall flutter (or near stall 
flutter). For the angle of attack and Mach 
number range considered (0<a<15°. 
0.6 < 0.95), the aerodynamic damping, 

although low at some test conditions remains 
positive (stabilizing) with the possible excep- 
tion of M = 0.90, a » 13°. Thus, stall flutter 
is not indicated except possibly in the vicinity 
of Af = 0.9, a = 13°. 


It is noted that the parameter bco 2 /a 
(a = speed of sound) had a value varying from 
0.50 to 0.52 for 0.6 <M 0 0 < 0.95. (A rough rule of thumb is that a value greater than about 0.5 for 
this parameter is sometimes indicative of a stall-flutter free wing (ref. 8).) 


CONCLUDING REMARKS 


Experimental results have been presented for the transonic flutter and buffet behavior for a 
low-aspect-ratio wing having an NACA OOXX-64 airfoil shape with a 9.3 percent thickness at the 
root and a 6.7 percent thickness at the tip. The model was tested at Mach numbers from 0.6 
to 0.95, at angles of attack from zero to 15°, and at various dynamic pressures. The following 
results are noted: 

1 . Single-mode type flutter instabilities occurred at near zero angle of attack. These instabil- 
ities were extremely sensitive to both Mach number and Reynolds number and appear to be 
associated with the onset of shock stall. The results suggest that delta wing configurations might 
also be susceptible to the same type of instability. 
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2. With one possible exception, stall flutter was not observed over the Mach number and 
angle-of-attack ranges investigated. 

3. Conditions of maximum buffet intensity occurred at angles of attack ranging from 10° 
to 14°. Quantitative results for the peak dynamic bending moments due to buffet are given in a 
form suitable for predicting full scale loads. 


Ames Research Center 

National Aeronautics and Space Administration 
Moffett Field, Calif. 94035, August 27, 1973 
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